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ABSTRACT 

The illumination pattern (or emissivity profile) of the accretion disc due to the reflec- 
tion of X-rays in AGN can be understood in terms of relativistic effects on the rays 
propagating from a source in a corona surrounding the central black hole, both on 
their trajectories and on the accretion disc itself. Theoretical emissivity profiles due to 
isotropic point sources as well as simple extended geometries are computed in general 
relativistic ray tracing simulations performed on graphics processing units (GPUs). 
Such simulations assuming only general relativity naturally explain the accretion disc 
emissivity profiles determined from relativistically broadened emission lines which fall 
off steeply (with power law indices of between 6 and 8) over the inner regions of the 
disc, then flattening off to almost a constant before tending to a constant power law 
of index 3 over the outer disc. Simulations for a variety of source locations, extents 
and geometries show how the emissivity profiles depend on these properties, and when 
combined with reverberation time lags allow the location and extent of the primary 
X-ray source to be constrained. Comparing the emissivity profile determined from the 
broadened iron K emission line in spectra of IH 0707-495 obtained in January 2008 to 
theoretical emissivity profiles and applying constraints from reverberation lags suggest 
that there exists an extended region of primary X-ray emission located as low as 2 rg 
above the accretion disc, extending outwards to a radius of around 30 rg. 

Key words: accretion discs - black hole physics - line: profiles - X-rays: general. 
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1 INTRODUCTION 

When the accretion disc surrounding the central black hole 
' of an AGN (or for that matter a galactic black hole) is illumi- 
nated by a non-thermal X-ray source in a corona of hot par- 
ticles surrounding the black hole, in addition to the contin- 
uum spectrum from this source, a reflectio n spectrum arising 
from the accretion disc will be observed ([George FabianI 
Il99lh . The reflection spectrum will con tain a number o f 
emission lines and absorption features (|Ross et al.l Il999h . 
broadened by general relativistic effects in the proximity of 
the black hole and relativistic Doppler shifting due to the 
moti on of the ref lecting material in its orbit in the accretion 
disc ([Fabian et a l. 1989; Laor 1991), the most prominent of 
which is the Ka emission line o f iron at 6.4 keV in the rest 
frame of the emitting material ([Matt et al.lll997l V 

In order to understand the emergent reflection spec- 
trum, it is necessary to understand the illumination pattern 
of the accretion disc, that is its emissivity profile^ the re- 
flected power per unit are a as a function of location on the 
disc. Following iLaorl ([ 199 ill ), the emissivity, e(r), is defined as 
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the radial weighting of the reflected radiation from the ac- 
cretion disc, measured in the rest frame of the disc material, 
such that the observed spectrum is given by 



Fo{vo) = / e(re)/r ( — J T(re,g)dgredre 



(1) 



where Ir is the rest-frame reflection emerging from the accre- 
tion disc, folded through the transfer funtion, T(re, g) which 
projects the rays around the black hole to the observer, inte- 
grating over all redshifts g that will shift the emitted photons 
into the observed energy band. 

In the simplest case of a point source in flat, Euclidean 
spacetime, the emissivity at a point on the disc is propor- 
tional to the inverse-square of the distance from the source, 
multiplied by the cosine of the angle at which the ray hits 
the disc from the normal, giving a form at large radius 
out from the source. In general relativity, however, rays will 
be focussed towards the black hole (and the inner disc) so 
naively one would expect a steeper fall-off in emissivity with 
distance from the black hole. 

Initial models for the spectru m of re l ativis tically broad- 
ened emission lines such as that of iLaorl ([l99lh assume that 
the emissivity profile of the accretion disc follows a power 
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law of constant index. The classical case would be an index 
of 3, though the best fits to observed spectra may require 
a steeper index. iDovciak et aD (|2004l l present an extended 
framework for fitting the spectral profiles of relativistically 
broadened emission lines with a suite of models including 
KYLINE which can either assume an emissivity profile taking 
form of a power law or take a tabulated radial emissivity 
profile for the accretion disc. 

To reproduce the enhancement in emissivity on the 
inner regions of the disc due to relativist ic effects, in or- 
der to correctly model the observed ir on Ka emissi o n line 
in the Seyfert 1 galaxy MCG-6-30-15, iFabian et al.1 (|2002l l 
found that a once broken power law form was required for 
the emissivity profile in the laor2 broadened emission line 
model, with the emissivity profile falling off more steeply 
over the inner regions of the disc before tending to approx- 
imately the classically-expected index of around 3 over the 
outer p art of the disc. This approach was also taken by 
brenneman Revnoldsl (|2006l ) using the observed profile of 
the iron Ka spectral line to measure the spin of the black 
hole (the KERRDISK model) and in the relline model of 
iDauser et all ^Old). 

While the emissivity profile is typically approximated 
as a power law when fitting models to the reflection spec- 
trum, the details will depend upon the location, spatial ex- 
tent and geometry of the source. A number of properties of 
accreting black holes are determined through the detection 
of emission from innermost regions of the accretion disc, 
not least the innermost extent of the disc which is used 
to infer the spin of the central black hole, assuming the 
accretion disc extends do wn to the innermost stable orbit 
(|Revnolds FabianI [iooi l or to infer truncation of the ac- 
cretion disc before this point. It is therefore important to 
understand how the emissivity profile of the disc varies over 
this region such that detections of emission from these radii 
can be correctly interpreted. For example, if too steep a 
falling emissivity profile is assumed, in order to correctly re- 
produce the flux levels emitted from the innermost regions, 
the disc must be truncated at a larger radius, in turn causing 
the spin of the black hole to be underestimated. 

Previous calculations of the illumination of black 
hole accretion discs by a poi nt source on the rot ation 
axis above the b l ack h ole, e.g. iMiniutti et al. 1 (|2003 h and 
iMartocchia et al.l (|2000l l, predict emissivity profiles approx- 
imating twice-broken power laws, falling steeply over the 
inner regions of the disc, then flattening before tending to a 
constant power la w inde x over the outermost radii, while 
ISuebsuwong et ahl (|2006l ) extended this work to orbiting 
point sources offset from the rotation axis, yielding quali- 
tatively similar results. 

The emissivity profile of the accretion disc in the nar- 
row line Seyfert 1 galaxy, IH 0707-495 was determined by in- 
dependently fitting the contributions to the relativistically 
broadened emission lines from successive radii in the disc 
(fwTlkins Fabianll201lh . The accretion disc was found to 
have a steeply falling emissivity profile (with a power law 
index 7.8) in the inner region of the disc then fiattening off 
to almost constant emissivity between radii of 5.6 rg and 
34.5 Tg before tending to a constant power law index of 3.3 
(slightly steeper than the classical case) over the outermost 
parts of the disc, in broad agreement with previously pre- 
dicted forms. 



We develop here a formalism in which theoretical emis- 
sivity profiles due to point sources and extended sources 
located in the corona may be predicted through general rel- 
ativist ic ray tracing simulations. While previous theoretical 
work derives emissivity profiles of accretion discs for selected 
cases, mostly for the case of a point sources on the rotation 
axis above the central black hole, we systematically calcu- 
late emissivity profiles for varying locations of the primary 
X-ray source and further this, to include extended sources 
of different sizes and geometries. This allows the effects of 
these parameters on the emissivity profile to be studied, such 
that properties of the X-ray source may be inferred from the 
observed emissivity profile. 



2 CALCULATING THEORETICAL 
EMISSIVITY PROFILES 

Theoretical emissivity profiles due to the refiection of the 
primary X-ray continuum, from a source in a corona in the 
central regions of the AGN, off of an accretion disc in the 
equatorial plane are calculated in ray tracing simulations in 
the Kerr spacetime around the black hole. 



2.1 Ray tracing 

In B oyer-Lindquist co-ordinates, the Kerr metric (i Kerd 
Il963h for the spacetime around a spinning black hole (spin 
parameter a = is written 

/ 2 , 2 , 2a^rsin^6'\ .2^,2 
— I r + a H 1 sm OdLp 

where 

2 2 I 2 2/1 

p = r -\- a cos u 
A - 2r + 

Taking c = 1 and = 1 to work in the natural units of 
gravitational radii, Vg — for both spacelike and timelike 
co-ordinates. 

In this spacetime, photons propagate along null 
geodesies satisfying the condition 

g^ux^x"" = (2) 

and the geodesic equation 

+ rVi^^i^' = (3) 

where gij,u is the metric tensor defined such that ds^ = 
g^jjX^x^ ^ and F^^^^ ^ire the connection co-efficients, defining 
the derivatives of the basis vectors with = F^^^^ec. Re- 
peated upper and and lower indices imply summation. 

Solving Equation [3] with constraint [21 the geodesic equa- 
tions that characterise the propagation of photons can be 
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written r = (t,r,0, cp) 

[{r^ + cos^ 0){r^ + a^) + 2a^rsm'^ O] k - 2arh 
r2 f 1 + - 2 

2arA: sin^ 6* + (r^ + cos^ - 2r)h 
(r2 + a2) (r2 + cos^ 6* - 2r) sin^ 6* + 2a^r sin^ 6' 
Q + (/ca cos ^ — /i cot ^) {ka cos ^ + /i cot ^) 



t 



LP 



6>2 



The derivatives are with respect to some affine parameter, a 
and A;, h and Q are conser ved con stants of the motion which 
distinguish the geodesies (|Carterl ir968). 

Given the starting point of a photon and the values of 
the constants of motion, a ray can be traced as the affine 
parameter advances by integration of the geodesic equations. 



r(cr + da) — r(cr) + r{a)da 



(4) 



A number of singularities exist in the Kerr spacetime 
when expressed in Boyer-Lindquist co-ordinates. In addi- 
tion to the co-ordinate singularity defining the event hori- 
zon, from which rays cannot return once they have crossed at 
rn = I — y^(l — a)(l + a) for spin parameter a, further co- 
ordinate singularities exist at ^ = and at cp — 0,27v (where 
the co-ordinate is discontinuous) . It is necessary to integrate 
the geodesic equation with greater numerical precision closer 
to the singularities, to prevent the solution diverging here, 
but for economy in the computation, a larger step size can 
be used further from these singularities. As such, a variable 
step size in the affine parameter is used, proportional to the 
distance from the singularity divided by the velocity in that 
co-ordinate, given a baseline precision parameter r, defined 
such that 



da — 



(5) 



The appropriate step size is computed in this manner for 
each of the r, and Lp co-ordinates, and the step size is 
taken to be whichever of these is the smallest. 



2.2 The X-ray source 

An isotropic point source is defined as emitting equal power 
into equal solid angle in its own instantaneous rest frame, 
where the solid angle element for polar angles a and f5 is 

dO! — d(cosa)dl3 

The source frame is defined by a tetrad of orthogonal basis 
vectors. The timelike tetrad basis vector, e^^^ is, by definition 
of the instantaneous rest frame, parallel to the source's 4- 
velocity, uo (since in its rest frame, the 4- velocity must have 
no spacelike component as its 3- velocity is zero). By the 
equivalence principle which states that in a local, freely- 
falling laboratory, the laws of physics reduce to those of 
special relativity, the local spacetime in the source frame is 
flat and the tetrad (unit) basis vectors satisfy 



(6) 



Note that bracketed subscripts indicate directions in the 
tetrad basis and those not bracketed indicate directions in 



the co-ordinate basis. ?7(a)(6) = diag(l, — 1, — 1, — 1) is the 
Minkowski (fiat) space metric. The 2;-axis of the source 
frame is constructed parallel to the radial basis vector in 
the (Boyer-Lindquist) co-ordinate frame (see Appendix). 

In the source frame, rays with photon energy E are 
started at equal intervals in solid angle, that is equal inter- 
vals in cos a and /3, with 4-momentum 



p = {E, E sin a sin /3, E sin a cos /3, E cos a) 



(7) 



The ray is then transformed back to the global co- 
ordinate basis and given the source location, the constants 
of motion for the ray can be found by inverting the equations 
of motion, allowing the ray to be propagated. 

2.3 Emissivity Profile Calculation 

Photons are propagated until they reach the disc, when ^ 
^ , or until they escape to a maximum allowed radius (or the 
maximum number of allowed steps is reached, to prevent the 
code entering an infinite loop if the step size were to become 
very small). When the photon hits the disc, its position is 
assigned to a radial bin and the number of photons hitting 
the disc in each bin is counted. 

Once the number of photons in each radial bin has been 
counted, the emissivity profile is obtained by dividing these 
counts by the area of each bin (an annulus). Strictly speak- 
ing, that is the proper area of the annulus, as measured in 
the disc frame, from the definition of the emissivity profile. 
The proper area as measured by a stationary observer is ob- 
tained from the metric. For radial bins of the accretion disc, 
dt — dO — 0^ so the area element is 

d^x — ^grrQipipdrdcp 



^ / + + ^^^^drdip 



(8) 



So the area of an annulus of co-ordinate thickness dr at 
radial co-ordinate r is 



A{r,dr) ^tt-^ 



2,2, 2a2r 
+ H —ar 



(9) 



In addition to the warping of spacetime around the 
black hole, the proper area of an orbiting element of the 
disc will be Lorentz contracted according to a stationary 
observer, so in the disc frame, the area will be increased by 
a factor of 7, the Lorentz factor of the orbiting disc element. 

A further subtlety lies here, and that is the velocity 
of the 'stationary' observer mentioned above. In the Kerr 
spacetime, the rotation of the black hole causes 'frame- 
dragging' of an observer whose angular momentum is zero. 
Therefore it is more meaningful to define the stationary 
observer as one who is dragged in this way, rotating at 
^ = = ^ (where = (r' + a')' - a^Asin^O). In 
fact, the solution ^ = is not possible within the station- 
ary limit surface upon which gu = (this lies at 2 rg for a 
maximally rotating Kerr black hole) — massive particles are 
forced to orbit the black hole within this region since there 
is no allowed solution of the geodesic equations which have 
zero 3- velocity. 

The Lorentz factor of the orbiting disc element is cal- 
culated from the disc velocity with respect to an observer 
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with = cj. Solving the geodesic equations for a massive 
particle in a circular orbit (r = 0) in the equatorial plane, 
the 4- velocity of the disc element is 

vd = (t,0,0,(^) = t(l,0,0,Q) 

(with the + and — signs in the denominator correspond- 
ing to prograde and retrograde orbits respectively). This is 
then projected into the frame of the observer (whose tetrad 
basis is constructed as for the X-ray source) by taking the 
appropriate scalar products. 

(a) I 11 V 

The Lorentz factor is then calculated from the 3-velocity, 
~^ — (v'^^^ ^v^^^ ^v^^^) and its squared magnitude, — 

(„(l))2 + („(2))2 + („(3))2. 

In addition to the relativistic effects on the areas of 
the radial bins, the energy of individual rays will be red or 
blueshifted as rays travel further from or closer to the black 
hole due to the variation in the rates at which the proper 
times of observers elapse. Since the emissivity is defined as 
the flux emitted from the disc (proportional to the flux re- 
ceived for reflection) and the flux is the power per unit area, 
it is proportional to the product of the photon arrival rate 
and the energy of each photon. The emissivity is enhanced 
by a factor oi ^ where ^ = is the ratio of the emitted 
and observed photon energies. The first factor of g arises 
from the shifts in the energy of individual photons, while 
the second is due to the photon arrival rate along each ray. 
The redshift is calculated by the projection of the photon 
4-momentum on its geodesic on to the observers' timelike 
axes since the photon energy is the p° component of the 
4-momentum. 

-1 ^ ^ ^ VQ • P(Q) ^ gy^v^oV'^iP) 

VE-p(E) gp.v'j,p-{E) ^ ^ 

where vq and ve are, respectively, the 4- velocities of the 
observer (the disc element) and emitter (the primary X-ray 
source), while p(x) is the photon 4-momentum as a function 
of position along the geodesic, here taken at observation and 
emission. 

Putting this together, if the photon count in a radial 
bin at co-ordinate radius r and of thickness dr is N(r,dr), 
the emissivity profile is given by 

^ N{r,dr) 

g^A{r,dr) ^ ^ 

2.4 Implementation 

The above algorithm is coded to run on Graphics Processing 
Units (GPUs) using the NVIDIA CUDA programming architec- 
ture. A typical GPU at the time of writing has between 256 
and 512 processing cores, each capable of running an inde- 
pendent computing thread, meaning that parallelised com- 
putations can often be sped up by factors of several hundred 
or more. In the CUDA programming model, threads are con- 
ceptualised on a two-dimensional grid, meaning each thread 
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Figure 1. Theoretical accretion disc emissivity profile due to a 
stationary isotropic point source located at a height of 10 rg on 
the rotation axis. 

can be readily mapped on to a point in the source-frame 
trajectory (cosa,/3) parameter space. 

In order to compute emissivity profiles due to an X- 
ray source from a given set of rays, each ray is traced by an 
independent computing thread on the GPU. The parameters 
(starting co-ordinates and constant of motion) are computed 
by the conventional CPU and set up in the memory on board 
the GPU, with memory addresses allocated to each thread 
{e.g. for an isotropic point source, all the rays are started 
at the same spacetime location and the constants of motion 
are incremented according to equal steps in cos a and /3). 
The computing threads are then executed, moving each ray 
as the affine parameter advances according to the equations 
of motion, until they reach their limits on the disc, at the 
event horizon or at a maximum allowed radius or number of 
steps. Once each thread has finished, the final positions of 
the rays are read back from the GPU memory and binned 
into locations on the accretion disc allowing the emissivity 
profile to computed. 



3 THEORETICAL EMISSIVITY PROFILES 
3.1 Relativistic Effects 

Computing the emissivity profile resulting from an isotropic 
point source at a height of 10 rg on the rotation axis above 
the plane of the accretion disc (Fig. [1]) illustrates the key 
effects that influence the form of the emissivity profile that 
will enable observed forms to be explained in terms of the 
properties of the X-ray source. 

In a flat, Euclidean spacetime the flux received from 
the source at each point on the disc and thus the reflected 
flux (emissivity) from that point will vary simply as the 
inverse square of the distance from the primary X-ray source, 
projected into the direction normal to the disc plane. If the 
source is at a height h above r = in the disc, the emissivity 
profile will go as (r^ + /i^)~"^ cosi^, with cos-i^ = , ^ , the 

Y r2-|-h,2 

angle from the normal at which the ray hits the disc. This 
will be constant (a flat profile) in the limit r <^ h (the inner 
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Classical area, r dr 
GR Disc Area 
Y * GR Disc Area 
r dr/ g 

y * GR Disc Area /g^ 



10 



100 



Figure 2. The effective area of successive annuli in the accretion 
disc, that is the total value (per dr) by which the photon count in 
the radial bin from the ray tracing simulation is divided in order 
to obtain the emissivity profile in Equation llll The classical area, 
rdr of the annulus is compared with that in the Kerr spacetime 
for a stationary area element as well as an orbiting element, mul- 
tiplying by the Lorentz factor, 7, to take into account Lorentz 
contraction as observed by a stationary observer. Finally, the ef- 
fect of redshift on the rays is accounted for. It is noted that the 
relativistic effects on the area of the disc element exactly cancel 
one factor of the redshift, gf, such that the fully relativistic result 
is equal to rdr j g. 



disc) and tending to at large radius, giving a power 
law emissivity profile with index 3 as can be seen where the 
flatter region of the emissivity profile becomes a power law of 
constant index at the break point, r ^ lOrg (approximately 
the height of the X-ray source above the accretion disc). 

The index of the power law over the outer region of 
the disc is, however, slightly steeper than the classical case 
of 3.0, with a value of around 3.2. In the presence of the 
black hole, gravitational light-bending will act to focus light 
rays from the point source towards the black hole and will 
therefore increase the flux incident on the inner regions of 
the disc relative to the outer regions causing it to fall off 
faster, steepening the emissivity profile. 

The contributions of the relativistic effects on the accre- 
tion disc itself are best illustrated by considering the 'effec- 
tive area' of the annuli in the accretion disc, that is the total 
value (per dr) by which the photon count in the radial bin 
from the ray tracing simulation is divided in order to obtain 
the emissivity profile in Equationllll Fig. [2] shows the classi- 
cal area, rdr of the annulus, compared with the proper area 
in the Kerr spacetime for a stationary area element as well 
as an orbiting element, multiplying by the Lorentz factor, 
7, to take into account Lorentz contraction as observed by 
a stationary observer. 

On including the shift in photon energy and arrival rate 
(factor of from Equation llOp as rays travel towards the 
black hole, the cause of the steep emissivity profile over the 
central regions of the disc is apparent. The increase in pho- 
ton energy as well as the increased photon arrival rate along 
each ray as measured by observers closer to the black hole 
with more slowly elapsing proper times, means that the flux 




Figure 3. The geometry of an isotropic point source located on 
the rotation axis above the black hole for which theoretical emis- 
sivity profiles are initially computed. 



received by the disc is greatly enhanced on the inner parts, 
increasing reflection from these regions (notwithstanding the 
effects of radiation transport out to the observer at infinity, 
as the emissivity profile is defined in the local rest frame of 
the emitting material before the transfer to the observer is 
included in the formation of the relativistic reflection spec- 
trum). 

It is interesting to note that the the general relativistic 
effects on the area of orbiting annuli in the accretion disc are 
cancelled exactly by one factor of redshift, so in practice the 
effective area of the annulus is given by dividing the classical 
area (rdr) by one factor of the redshift. This function is also 
plotted for comparison. 



3.2 Axial X-ray Sources 

The simplest, idealised case is that of an isotropic point 
source, stationary upon the rotation axis (Fig. [3]). This will 
allow us to explore the effects of ray propagation in the 
Kerr spacetime on the emissivity profile of the accretion 
disc, with the fewest free parameters and assumptions about 
the nature of the source itself. A localised source, however, 
may be expected if the X-ray emission re sults from magnetic 
reconnection events within the corona (|Galeev et gllll979l : 
iMerloni h Fabianll2QQlh . 

Theoretical emissivity profiles for isotropic point 
sources, stationary upon the rotation axis, at varying heights 
above the black hole are shown in Fig. [H 

In all cases, the profile tends to a power law with an 
index slightly steeper than 3 over the outer regions of the 
disc (from 3.1 for a source at height 10 rg to a steeper index 
of 3.3 for a source height of 3 rg), with the power law indices 
steepening as high as 6 ~ 7 over the innermost parts. 

It can be seen that as the source is moved higher up 
the rotation axis, further from the black hole, the region 
over which the emissivity profile is flattened increases (the 
region where r <^ h). The outer break-point in the power law 
form of the profile moves to coincide approximately in radius 
with the height of the source (Fig. [5|) for sources at heights 
greater than 12 rg with the outer break at a radius slightly 
greater than the height of the source for lower sources (due 
to the requirement for the break being r ^ h which requires 
a greater increment in r for smaller h). 

As the source is moved closer to the black hole, the 
steepening over the inner regions is greatly enhanced as more 
photons are focussed on to the inner regions of the disc while 
rays are bent towards the central black hole, until the source 
is as low as 3 rg at which point the steepening is so significant 
that it masks the classically-predicted flattened region of the 
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Figure 4. | (a) | Theoretical accretion disc emissivity profiles due to a stationary isotropic point sources located at varying heights 
above the black hole on the rotation axis. |(b)| The emissivity profiles scaled by r^, illustrating the steepening of the emissivity 
profile over the inner disc (where the plot decreases) before the profile flattens (plot increasing) and then tends to a constant power 
law index of around 3 over the outer disc (plot approximately constant). The plot styles in the two plots correspond to one another. 
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Figure 5. Observed location of the outer break point between 
the flattened region of the emissivity profile and the power law of 
constant index around 3 over the outer part of the accretion disc 
(where r ^ h) for axial point sources of increasing height above 
the disc. 



profile (now constrained to a much smaller region of the disc 
where r <^h). The emissivity profile more closely resembles 
a once-broken power law which in the extreme case of a 
source at only 1.235 rg above the black hole has a very steep 
power law index of around 8 out as far as a radius of 5rg 
before tending towards an index around 3.5 over the outer 
regions of the disc. 



3.3 Orbiting Sources 

While it is conceivable that X-rays originate from a station- 
ary point source close to the rotation axis, one might also 
consider the case of a point source located elsewhere in the 
corona. The concept of a stationary object in the corona (if 




Figure 6. The considered geometry of an isotropic point source 
orbiting the rotation axis above the plane of the accretion disc. 



not on the rotation axis) is somewhat unphysical, given that 
the AGN is considered to have formed from the gravitational 
collapse of material towards the galactic centre, a process in 
which it is likely to rotate, conserving angular momentum. 
Furthermore, without rotation, the material will just fall 
into the black hole and will not survive for long as a corona 
(unless it is replenished). Due to the axisymmetry of the 
Kerr spacetime, a point source in orbit at a given radius is 
equivalent to a continuous ring source of that radius (which 
may be relevant when considering the total emission from 
many localised flaring events when looking at time- averaged 
X-ray spectra). 

We again consider the idealised case of an isotropic 
point source to explore the observed effects in the emissivity 
profile due to ray propagation to the accretion disc in gen- 
eral relativity with the minimum number of free parameters. 
Orbiting point sources at various locations in the corona will 
be the building blocks for extended X-ray sources. 

Theoretical emissivity profiles for isotropic point 
sources as shown in Fig. [6l orbiting the rotation axis at 
varying radii at a height h — 5rg above the disc plane are 
shown in Fig. [71 The sources are 'co-rotating' with the ele- 
ment of the accretion disc at the same radius, i.e. a source 
at a distance x from the rotation axis as measured along 
a plane parallel to the accretion disc below is taken to be 
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Figure 7. Theoretical emissivity profiles for isotropic point 
sources orbiting the rotation axis at varying radii at a height 
h = 5rg above the disc plane. The sources are 'co-rotating' with 
the element of the accretion disc at the same radius. 

orbiting at the same velocity as the element of the disc in a 
(relativistic) Keplerian orbit at radius x. Where the orbiting 
source is close to the disc, this serves as an approximation to 
the orbital velocity, however if the X-ray source originates 
from flaring due to magnetic reconnection in poloidal field 
lines anchored to the ionised accretion disc, it may also be 
expected that the coronal material will move along with the 
orbiting disc. 

As for the case of an axial source, the profiles are steep- 
ened from the classical case due to gravitational light bend- 
ing, focussing more rays on to the inner region of the disc as 
well as the blueshifting of photons as they travel towards the 
innermost regions. These calculations show that the power 
law index of the emissivity profile over the innermost regions 
of the disc ranges from 7 for sources close to the rotation axis 
and black hole, within 5rg, to 6 where the source located 
further out between 20 and 25 rg (this central steepening of 
the emissivity profile is almost entirely due to the blueshift- 
ing of photon energies and arrival rates as the rays travel 
towards the black hole now the source is located away from 
the inner parts of the disc) . The index of the emissivity pro- 
file for the outer disc is 3.2 in each case. 

The emissivity profiles are flatter over the middle region 
than for axial sources, owing to the source being located 
further out from the rotation axis, so the classical region 
where h ^ r giving a constant flux at the disc is not masked 
by the relativistic effects steepening the profile over the inner 
disc. 

Since an X-ray source close to the innermost stable or- 
bit will be travelling at relativistic speeds if in a Keplerian 
orbit, emission will be 'beamed' into the forward direction 
of source motion. This will enhance the flux in front of the 
source and (to a lesser extent) on the region of the disc di- 
rectly below the source locus while reducing the flux received 
outside of this region. This will serve to further flatten the 
emissivity profile. 

Following this reasoning, increasing the radius of the 
ring source flattens the profile out to a larger radius. Where 
the source height is less than or of the order of the source's 



radius, a small peak is observed in the emissivity profile 
below the source, where the proximate regions of the disc 
subtend a large solid angle at the source so a large number 
of rays are intercepted by this region of the disc before they 
are able to propagate further out. 

3.4 Extended Sources 

While it is instructive to consider the case of isotropic point 
sources in the corona to illustrate the factors affecting the 
emissivity profile of the accretion disc from X-ray reflection, 
in reality, the source is likely to be extended over a region 
of the corona where particles are heated and accelerated, 
up-scattering seed photons to the observed X-rays. 

In the simplest case, an extended X-ray source can be 
considered as the sum of point sources assuming the corona 
is optically thin to the emitted X-rays such that emission 
from any part of the source can escape without interacting 
with other parts (as one might expect for an optically thin 
corona of hot, accelerated electrons that inverse- Compt on 
scatter photons). This simple model will illustrate the be- 
haviour of such systems in determining the accretion disc 
emissivity profile. 

A vertically extended source of X-ray emission can be 
constructed by summing isotropic point sources at regular 
intervals. This will allow the emissivity profiles of sources 
of finite height above the disc plane to be investigated. The 
emissivity profile resulting from a stationary, vertically ex- 
tended source whose luminosity is constant along its length 
is shown in Fig. (8] These profiles take an almost constant 
power law index all the way out on the disc as the profiles 
from sources at successive heights sum together, masking 
the flattened part of each. There is steepening due to rel- 
ativistic blueshift on the innermost parts and a very slight 
break to the outermost index of around 3 corresponding to 
the break point in the emissivity profile of the uppermost 
extent of the source. 

Spatially extended X-ray sources may be studied in 
more detail through Monte- Carlo ray tracing simulations. 
Rays are started at random locations (with a uniform prob- 
ability density function such that all locations are equally 
likely) within an allowed cylindrical source region (defined 
by a lower and upper height from the disc plane as well as 
an inner and outer radius) and are assigned random initial 
direction cosines cos a and /3, again with uniform probabil- 
ity density (Fig. [9l This will simulate the effect of an X-ray 
source of finite spatial extent that is, again, optically thin 
to the X-rays it emits. Each local region of the source is 
taken to be co-rotating with the element of the accretion 
disc in a relativistic Keplerian orbit at the same radius. To 
account for the variation in source luminosity across its ex- 
tent, the rays reaching the disc can be weighted by power 
laws in the height and radius of their origin in the source (or 
indeed weighted by other functions) allowing for simple non- 
uniform sources to be modelled, however as this introduces 
a number of free parameters that will not be constrained by 
the current quality of observations, the simplest case of con- 
stant luminosity sources is taken here to explore the effects 
of spatially extending the source. 

Modelling a radially extended source at a height of 10 rg 
and extending radially to 25 rg produces an emissivity profile 
as shown in Fig. [10] and, for reference, is compared to a 
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Figure 8. Theoretical accretion disc emissivity profile due to a 
vertically extended stationary source on the rotation axis, extend- 
ing from 2 Tg to 20 rg. The emissivity profile compared to those of 
point sources located at either end of the extended source, illus- 
trating that the overall emissivity profile combines the effects of 
the more extreme steepening over the inner disc from the source 
closer to the black hole, with the slight outer break point cor- 
responding to that in the emissivity profile of the highest source 
(though the effect of summing the sources up the axis is to steepen 
the middle region of the profile, making the break points in the 
power law form less pronounced). 




Figure 9. Extended X-ray source defined by a lower and upper 
height from the disc plane as well as an inner and outer radius. 



single point source located at the outer extent. It can be 
seen that the outermost break point in the power law form 
is determined by the outer extent of the source, while the 
existence of X-ray sources within this radius cause the profile 
to be flattened off within this part before steepening, once 
again, over the inner regions of the disc. 



3.5 Jet Sources 

X-ray emitting jets of collimated particles accelerated to rel- 
ativistic velocities are observed in a number of galactic black 
hole binaries such as Cygnus X-1. It is conceivable that it 
is the X-rays from this jet that illuminate the accretion disc 
giving rise to the reflection dominated component of the 
spectrum. A vertically collimated X-ray source could also 
represent, for example, the X-ray emission from a vertically 
collimated jet of particles accelerated up along the rotation 
axis of the black hole in a radio galaxy. 



Figure 10. Theoretical accretion disc emissivity profiles due to 
an extended disc of emission located at a height of 10 rg above 
the disc plane and extending radially to 25 rg from the rotation 
axis, compared to that arising from a point source orbiting on the 
outermost edge of this source. 



dr/dt = 0(V=0) 
dr/dt = 0.1c (V= 0.16 c) 
dr/dt = 0.2c (V= 0.33 c) 
dr/dt = 0.3c (V= 0.49 c) 
dr/dt = 0.4c (V=0.65 c) 
dr/dt = 0.5c (V= 0.81 c) 




r/rg 



Figure 11. Theoretical accretion disc emissivity profiles due to a 
'jet' component moving radially up the rotation axis. The moving 
point source is taken to be at a height of 5 rg and is assumed to be 
isotropic in its instantaneous rest frame. Velocities are quoted in 
the global Boyer-Lindquist co-ordinates (dr/dt) as well as those 
measured by a stationary observer at the same location (V). The 
rays propagating backwards, towards the black hole, are still fo- 
cussed onto the inner parts of the accretion disc, however as the 
relativistic motion of the source outwards causes emission to be 
beamed away from the accretion disc causing a sudden drop in 
the emission reaching the disc over the middle region. 



Fig. \TT\ shows the accretion disc emissivity profiles due 
to illumination by point sources moving radially at a con- 
stant velocity up the rotation axis (the point sources are 
assumed to be isotropic in their instantaneous rest frames) 
compared to a stationary source. The emissivity profile are 
computed, as before, by constructing the tetrad basis vec- 
tors in the instantaneous rest frame of the X-ray source, now 
moving radially. 
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Figure 12. Relativistic motion of the jet source causes the emis- 
sion to be beamed into the forward direction of motion. The black 
hole still focusses the backward emission onto the innermost re- 
gions of the accretion disc, however X-rays are now beamed away 
from the accretion disc, reducing reflection from the middle re- 
gion. 

Gravitational light bending and blueshifting of rays still 
enhances the emission reaching the very inner part of the 
disc, however the highly relativistic radial motion of the 
sources causes the emission to be beamed in front of the 
motion, greatly reducing the emission behind the jet that 
reaches the middle region of the accretion disc (Fig. I12|) . 
The emission reaching the accretion disc drops dramatically 
by 1.5 ~ 2 orders of magnitude on the middle parts. 

The very steep fall off in the emissivity profile, the form 
of the profile greatly differing from observed emissivity pro- 
files and the low fraction of the initial radiation that is re- 
flected suggests that relativist ically moving jet sources are 
not responsible for the signiflcant reflection components ob- 
served in AGN and galactic black hole binaries, rather these 
more likely occur from a more slowly moving corona. 

Accelerating jets of particles do, however, raise the pos- 
sibility that the observed X-ray continuum and reflected 
component could be disconnected, with a fairly constant re- 
flection component from a steady base of a jet or extended 
corona while the observed variable continuum radiation is 
dominated by fast moving particles in the jet accelerated 
by a variable mechanism. Weaker variation of the reflec- 
tion component that is uncorrelated with the variation in 
the pr imary continuum was obser v ed bylVaughan EdelsonI 
(|200lh and iFabian k VaughanI (|2003l ) in MCG-6-30-15. 
Such disconnection of the observed continuum and reflection 
also relaxes the constraint on the relative fluxes observed in 
the continuum and reflection components of the spectrum if 
they are no longer dominated by the same component. 

3.6 Black Hole Spin 

The most significant effect of varying the spin of the black 
hole is to change the location of the innermost stable circu- 
lar orbit (ISCO). For a maximally rotating black hole (a = 
0.998) in the Kerr spacetime, the ISCO lies at a radial co- 
ordinate of 1.235 Tg while for a non-rotating, Schwarzschild 
black hole, the ISCO moves out to 6rg. The accretion disc 
cannot exist stably within the ISCO since material is unable 
to maintain a stable orbit here and will plunge into the black 
hole, reducing the density in this region causing little or no 
refiection to be seen within the ISCO. 

Fig. [13] shows theoretical emissivity profiles for the ac- 
cretion disc around black holes with varying (dimension- 
less) spin parameter, a, for both an axial and orbiting ring 
sources. It can be seen that the spin parameter has little 
effect on the emissivity at specific locations upon the disc. 



with the only notable change in the profiles being the trun- 
cation of the disc at greater radii as the spin parameter 
decreases. As such, the steepening over the inner part of the 
accretion disc is only seen for rapidly spinning black holes, 
with no steepened inner part seen for a ^ 0.8, though this is 
due to the lack of refiector at small radius where the spin is 
low rather than an intrinsic effect of the black hole spin. In 
the case of the source orbiting with the same angular veloc- 
ity as the disc element below, refiection is slightly enhanced 
from the innermost regions for lower values of the spin owing 
to the angular velocity of a circular orbit being greater for 
smaller values of the spin parameter, increasing relativistic 
beaming. 

This is, of course, only accounting for the propagation 
of X-rays around the black hole and assuming 'idealised' re- 
fiection off of a razor thin accretion disc extending as close 
to the black hole as it is able to as defined by the innermost 
stable orbit. It may be that the accretion disc does not ex- 
tend this far inwards, in the case of a truncated accretion 
disc with an advection-dominated hot accretion fiow in the 
inner regions, though this solution would be more applicable 
to accretion discs around galactic black holes rather than to 
AGN. For example, Beckwith et al. (2008) find in general 
relativistic magneto-hydrodynamic simulations that there is 
significant thermal dissipation within the innermost stable 
orbit, regardless of the spin of the black hole. The density 
of the accretion fiow will drop greatly in this region (both 
due to the high temperatures and there being no stable cir- 
cular orbit such that material plunges into the black hole), 
so little refiection, with which we are concerned here, will 
be observed from the accretion fiow within the innermost 
stable orbit. 



3.7 Returning Radiation 

In general relativity, it is conceivable that X-rays refiected 
from one part of the accretion disc will be bent around the 
black hole back to t he disc to be re fiected a second time 
(Cunningham 1976). Ross et al. (2002) considered the spec- 
trum of the refiected radiation from an ionised slab [i.e. an 
accretion disc) if a fraction of the refiected radiation is said 
to return to the refiector as incident radiation. They found 
that multiple refiection of returning radiation strengthens 
the broadened absorption and emission features such as the 
iron K line as well as steepening the softer part of the spec- 
trum below 2keV. We here simply consider the effects of 
gravitational light bending returning radiation to the disc 
in terms of the 'raw' emissivity rather than considering the 
energy-dependent details of the refiection spectrum, since 
the emissivity profile is determined from the emission lines 
which are enhanced so we determine if this measurement is 
affected by returning radiation. 

In order to compute the emissivity profile due to radi- 
ation returning to the accretion disc from the primary re- 
fiection component, the emission from the disc is modelled 
as the sum of isotropic point sources in (relativistic) Kep- 
lerian orbits in the disc at regular radial intervals. Due to 
the axisymmetry of the Kerr spacetime, each of these point 
sources represents emission (or more specifically refiection) 
from its respective annulus. The photon count from each 
source is scaled by the emissivity profile across the disc to 
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Figure 13. Theoretical emissivity profiles of the accretion disc surrounding black holes with varying (dimensionless) spin parameter, 
a for |(a)| a stationary source located at a height 5 rg on the rotation axis and |(b)| a ring source radius 3 rg located at a height 5 rg 
co-rotating with the disc below. 



account for the variation in incident flux and is multiplied 
by the (relativistic) area of the annulus. 

This, of course, assumes that in the frame of the disc 
material, the reflected radiation is emitted isotropically 
while in reality anisotropic emission from atomic processes 
and Compton scattering will dep end on the an ^le of inci- 
dence of the incoming radiation (Sv oboda et al.l l2010). how- 
ever provides a simple model which illustrates an upper limit 
of the effect of returning radiation. 

Simulations reveal that as many as 50 per cent of the 
photons reflected from the accretion disc are returned to the 
disc (since for a typical X-ray source located above the plane 
of the accretion disc, the reflected flux will be concentrated 
on the innermost parts of the disc where the influence of the 
black hole focusing them back to the disc is the greatest). 

Computing the returning radiation from a disc illumi- 
nated by a point source at a height 10 rg above the disc 
around a maximally spinning black hole yields an emissiv- 
ity profile with an approximately constant power law slope, 
steepening over the innermost regions of the disc. Once the 
photon counts returning to the disc have been divided by the 
effective areas of the annuli, the emissivity due to return- 
ing radiation is an order of magnitude lower than that due 
to the primary X-ray source and with its almost constant 
power law slope, the returning radiation has little effect on 
the overall shape of the emissivity profile, as shown in Fig. 

M 

These results are consistent with those o f lAgol Krdi^ 
(|2000l ) who show that for a maximally spinning black hole, 
up to 50 per cent of the radiation emitted from the accretion 
disc can be returned to the disc before finally being observed 
at infinity, though we here explicitly demonstrate that the 
distribution of re-refiected emission from the disc means 
that returning radiation does not grea tly affect the observed 
emissivity profile of the accretion disc. Agol & Krolik 
also note that the fraction of the flux returned to the accre- 
tion disc for more slowly spinning black holes is less (as many 
of the returning rays end up within the innermost stable or- 
bit in this case) so as returning radiation in the maximally 
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Figure 14. The effect of refiected radiation returning to the ac- 
cretion disc to be refiected again on the overall emissivity profile. 
The emissivity profiles due to the primary X-ray source and re- 
turning X-rays are shown as red crosses and green diamonds, re- 
spectively, and the overall profile (the sum of these contributions) 
is shown as black points. Returning radiation has little effect on 
the overall emissivity profile. 

spinning case does not affect the emissivity profile substan- 
tially, it will do so even less for more slowly spinning black 
holes. Furthermore, we can say that even including the ef- 
fects of returning radiation, the black hole spin does not 
greatly affect the form of the emissivity profile save for the 
innermost extent of the accretion disc and thus the degree 
of steepening observed over the inner part of the disc. 



4 THE CASE OF IH 0707-495 

The emissivity profile of the accretion disc due to X-ray re- 
flection in the narrow line Seyfert 1 gala xy, IH 0707-495 was 
determined bv I Wilkins &: FabianI (|201lh , considering the rel- 
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Figure 15. | (a) [ Emissivity profile determined for X-ray reflection 
from the accretion disc in IH 0707-495 in January 2008 by fitting 
for the relative contributions of components of the relativistically 
broadened iron K emiss ion line from successive radii in the disc 
(|Wilkins FabianllioTlh . compared with |(b)| a theoretical emis- 
sivity profile due to an extended X-ray source extending radially 
outwards to 30 rg and between 2 and 10 rg above the plane of the 
accretion disc. 

ativistically broadened iron K line (whose rest-frame en- 
ergy is 6.4 keV) to be composed of the sum of independent 
relativistically-blurred emission line components from suc- 
cessive annuli in the accretion disc and fitting for the relative 
contributions of these. 

The accretion disc was found to have an emissivity pro- 
file approximated by a twice-broken power law, with a steep 
index of 7.8 over the innermost parts then flattening to al- 
most constant emissivity (index zero) between 5.6 rg and 
34.8 Tg before tending to a constant index 3.3 over the outer 
part s of the disc (Fig. [15]). 

[Fabian et al.1 (|2009|) and lZoghbi et al.1 (|2010h report re- 
verberation time lags of around 30 s where variability in the 
reflection dominated component of the spectrum is seen to 
respond to that in the power law continuum after a time 
lag due to the finite light travel time from the primary X- 
ray source to the reflecting accretion disc. If the mass of 
the black hole is t aken to be 3 x 10^ as quoted in the 
literature (see e.g. IZhou &: Wan3l2005l V this implies that a 



significant portion of the observed X-rays originate from a 
source located within 2 rg of the reflector. 

Comparing the observed emissivity profile to the theo- 
retical emissivity profiles computed for a variety of source 
locations and geometries while applying the constraint that 
the source is as close to the disc as 2rg implies that the 
source is located at a low height close to the disc plane while 
extending radially outwards to between 25 and 35 rg to ac- 
count for the location of the outermost break point in the 
emissivity profile. A theoretical emissivity profile for a disc 
of emission between 2 rg and 10 rg above the disc, extending 
radially outward to 30 rg with constant luminosity through- 
out its extent is shown in Fig. \lE\ where letting the source 
region extend to greater heights serves to flatten the middle 
section of the emissivity profile. To account for the 30 s re- 
verberation signal, it is likely that the source will be more 
luminous in the lower regions closer to the disc, however 
this is also accounted for in the fact that a larger fraction of 
the primary emission will reach the disc from lower sources 
than those further from the disc, as can be seen in Fig. [4] al- 
though a fully relativistic treatment of reflection time lags is 
required to quantify this fully and will be explored in future 
work. 

The only notable discrepancy between theoretical mod- 
els and the observed emissivity profile of IH 0707-495 is the 
extent of the inner region over which the profile is steepened, 
where fitting to the profile of the iron K emission line shows 
the power law emissivity profile to break from the steep in- 
ner index to the flat region around a radius of 5rg, while 
the theoretical models produce emissivity profiles falling off 
steeply over the inner around 3 rg before the power law index 
gradually reduces to the flat region. In the observed spectra 
from which the emissivity profile of IH 0707-495 was deter- 
mined, relatively little flux is in the redshifted wing of the 
emission line since this originates from the innermost regions 
of the accretion disc closest to the black hole which subtend 
the smallest solid angle at the observer (not least because 
many of the rays reflected form this part of the disc will fall 
into the black hole). As such, the emissivity determined for 
this region is an upper limit, since to a certain degree, the 
emission from the inner parts of the disc can be increased 
with little effect seen by the observer as the excess photons 
fall into the black hole and are not observed. From the con- 
fidence limits on these contributions to the emission line, 
however, it is clear that some steepening of the emissivity 
profile is required to explain the broadened line profile in 
agreement with the theoretical calculations. 

The spectral features contributed to the emission line 
by the inner regions of the disc out to 10 r g will appear i n the 
energ y range 3-5 keV (see. e.g.. Fig. 1 of Wilkins FabianI 
1201 ih over which the intrinsic energy resolution (FWHM) of 
the EPIC pn detector, used to obtain the spectra from which 
this emissivity profie was determined, is 150 eV. Spectral fea- 
tures distinguishing the emission from these inner annuli of 
the accretion disc will, to some degree, be smoothed out 
by the detector such that they will be harder to distinguish 
above the power law continuum. The curving of the emissiv- 
ity profile to the flattened middle section is therefore likely 
not fully resolved, rather we fit a broken power law to the 
observed emissivity profile. 

The results of these ray tracing calculations combined 
with the observed emissivity profile of IH 0707-495 and the 
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reflection time lags appear to suggest an emitting region 
surrounding the central black hole and covering the cen- 
tral parts of the accretion disc. This emitting region could 
consist a corona of hot electrons emitting X-rays through 
the inverse- Compt on scattering of (thermal) seed photons 
emitted from the accretion disc. If this corona were opti- 
cally thin to X-rays, as implied by the steep index of the 
power law continuum after Comptonisation of seed pho- 
tons, for example a sparse plasma or even a region with 
emission origina ting from flares where magnetic reconnec- 
tion takes place (iGaleev et al.lll979l : iMerloni FabianlliooH : 
iGovder Lasenbvll2004i r rays emitted from any part of it 
could reach the accretion disc, as modelled in ray tracing 
simulations here and the corona would not obscure the emis- 
sion from the regions of the disc below. 

In January 2011, IH 0707-495 dropped into a low flux 
state in which the observed spectrum was dominated by the 
reflection from the accretion disc, with little to no direct 
emission seen from the primary X-ray source (|Fabian et al.l 
I2OI2I I . The emissivity profile of the accretion disc was deter- 
mined as for the January 2008 observation and was found 
to fall off steeply with a power law index of about 8 out to 
a radius of 5rg before tending to a constant index slightly 
steeper than the classical value of 3 over the outer part of the 
disc fFig. I16p . Comparing this to the theoretical emissivity 
profiles suggests that the X-ray source had collapsed down to 
a compact source around the rotation axis located at a low 
height of between 1.5 and 2 rg above the disc plane (Fig. [15]). 
Ray tracing simulations simply counting the number of pho- 
tons that land on the accretion disc compared to the number 
able to escape to infinity (which will be seen as the power 
law continuum in the observed spectrum) demonstrate that 
such a compact source close to the black hole naturally ex- 
plains the drop in the power law continuum emission with 
the majority of photons being focussed on to the accretion 
disc (and reflected) or into the black hole itself. 



5 CONCLUSIONS 

Numerical ray tracing simulations of the propagation of X- 
rays from a source in a corona surrounding the central black 
hole to the accretion disc have yielded theoretical predictions 
for the emissivity profile of the accretion disc due to X-ray 
reflection. 

These calculations have demonstrated that the forms 
of the emissivity profiles observed through the profiles of 
relatively broadened emission lines from accretion discs in 
AGN (such as that from the narrow line Seyfert 1 galaxy 
IH 0707-495) arise naturally from the relativistic effects on 
the propagation of the rays around the black hole and on the 
accretion disc itself, with the simplest possible assumptions 
about the X-ray itself (assuming either an isotropic point 
source or an extended source of uniform intensity through- 
out its volume). 

By comparing observed emissivity profiles to those com- 
puted theoretically for different locations and geometries of 
the source along with consideration of the time lags between 
the continuum emission and reflection components in AGN, 
it is possible to constrain the location and extent of the pri- 
mary X-ray source. In IH 0707-495, the emissivity profile 
obtained from the data appears to suggest an X-ray source 
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Figure 16. | (a) | Emissivity profile determined for X-ray reflection 
from the accretion disc in IH 0707-495 in January 2011 by fitting 
for the relative contributions of components of the relativistically 
broadened iron K e mission line from successive radii in the disc 
jFabian et al.ll20l3 ). compared with |(b)] a theoretical emissivity 
profile due to a compact X-ray source located on the X-ray source 
1.5 rg above the plane of the accretion disc. 



as low as 2rg above the plane of the accretion disc and ex- 
tending outwards from the rotation axis to around 30 rg. 

Considering the emissivity profile of the accretion disc 
in terms of the properties of the primary X-ray source has 
also allowed the change in state of 1H0707-495 observed in 
January 2011 (in which the spectrum was almost entirely 
reflection dominated with little direct continuum emission) 
to be interpreted in terms of the X-ray source collapsing to 
a confined region around the rotation axis, around 1.5-2 rg 
above the accretion disc. 
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Where 
e - 

With 



sin^ 



P_ 
A 

2iiacr 



S — (r + a 



2 \ 2 2 a •2/1 

- a A sm 6/ 



Hereafter, we shall set ji — c— 1. 

This definition of the metric co-efficients as exponentials 
allows their square roots, products and ratios to be easily 
written. 

The source frame is defined by a tetrad of orthogonal 
basis vectors. By the equivalence principle, which states that 
in a local, freely- falling laboratory the laws of physics re- 
duce to those of special relativity, the local spacetime in the 
source frame is flat and the tetrad basis vectors satisfy 



V{a){b) 
r](a){b) 



(Al) 
(A2) 



Note that bracketed subscripts indicate directions in the 
tetrad basis and those not bracketed directions in the 
co-ordinate basis. ?7(a)(b) = diag(l, — 1, — 1, — 1) is the 
Minkowski (flat) space metric. 

The timelike tetrad basis vector, e^^^^ is, by definition 
of the instantaneous rest frame, parallel to the source's 4- 
velocity, uo (since in its rest frame, the 4- velocity must have 
no spacelike component as its 3- velocity is zero). The 4- 
velocity of the emitting particle orbiting at angular velocity 



dt 



i^'(l,0,0,Q) 



The time- like component, ix*, is found from the condition for 
the motion of a massive particle g^^jU^u^ — c? (again taking 
c=l). 

gttuu + 2gti^uu^ + g^^u^u^ — 1 



2 2'0x/ t\2 . o 

- (jj e )[u ) + 2(jje 



2ip 



APPENDIX A: CONSTRUCTING THE 
SOURCE FRAME 

For completeness, we include here the derivation of the 
tetrad basis vectors defining the local rest frame of a source 
orbiting the rotation axis of the black hole in the Kerr space - 
time (see also, e.a.. iKrolik et al . 2005; Beck with et ahllioosh . 
When considering the dynamics of inertial frames and con- 
structing tetrad basis vectors, it is o ften convenient to re - 
write the metric in the form (following IChandrasekharll 1 9831 ) 

ds^ — e^^ c?dt^ — e^^ (dcp — uodt)^ — e^^^dr^ — e^^'^dO^ 



7l-e2(V^--)(ll-cj)2 

Since the 4-velocity has no component in either the r 
or directions and the Kerr metric is diagonal in these di- 
rections, two of the orthogonal space-like basis vectors are 

6(2) = (o,0,io) 

Finally, the fourth orthogonal tetrad basis vector must take 
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the form 

e(i) = (eji),0,0,e^,)) 
and satisfy the conditions ()A2p 

• 6(1) = (A3) 
6(1) -6(1) = -1 (A4) 
Condition ()A3|) gives 

^^6(^)6(1) +^tc^e(i)e^-L) +^c^te^t)e(i) + ^c^c^e^^^e^-^^ = 



^(1) - p2^p(o.-nA ""(1) 



And from (lA4l) . 

^"6(1)6(1) + 2^tc^e(i)e^-L) +^v'v'efi)efi) = -1 
Substituting using (|A|) and solving gives 



Hi) 



Where the signs have been selected to form a right-handed 
basis. 

Thus, the tetrad basis vectors in the frame of a massive 
observer orbiting a Kerr black hole at angular velocity ^ = 
Q are 



7l-e2(V^--)(Q-cj)2 



,0,0, 



^(1) 



7l-e2(V^--)(Q-cj)2 



,0,0, 



^g2. _g2^(Q_^)2 



3(2) = 



0,0,i,0 



-(3) = ^0,^,0,0 

Taking 6^3), 6^2) and e^^^ as the z (polar), x and y axes re- 
spectively (forming a right-handed co-ordinate system), the 
4-momentum of a photon energy Eq emitted at declination 
a and azimuth /S in the emitter's frame is written 

= {Eo , ^0 sin a sin /3, Eq sin a cos /3, Eq cos a) 

And transforming back into the Boyer-Lindquist co-ordinate 
basis, 

; (t) t , (1) t 

(3) r 

r = ^6(3) 

/3 (2) 

= p' ^6(2) 

(t) I (1) 

if = ^e^,) ^e^,) 
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